Determining the structure of Ru(0001) from 
low-energy electron diffraction of a single 

terrace 



J. de la Figuera, a J. M. Puerta, b J. I. Cerda, b F. El Gabaly, 

and K. F. McCarty c 

a Dpto. de Fisica de la Materia Condensada and Centro de Microandlisis de 
Materiales, Universidad Autonoma de Madrid, Madrid 28049, Spain 

h Instituto de Ciencia de Materiales de Madrid, CSIC, Madrid 28049, Spain 

c Sandia National Laboratories, Livermore, California 94551 



Abstract 

While a perfect hep (0001) surface has three-fold symmetry, the diffraction pat- 
terns commonly obtained are six-fold symmetric. This apparent change in symmetry 
occurs because on a stepped surface, the atomic layers on adjacent terraces are ro- 
tated by 180 degrees. Here we use a Low-Energy Electron Microscope to acquire 
the three-fold diffraction pattern from a single hep Ru terrace and measure the 
intensity-vs-energy curves for several diffracted beams. By means of multiple scat- 
tering calculations fitted to the experimental data with a Pendry R-factor of 0.077, 
we find that the surface is contracted by 3.5(=L0.9)% at 456 K. 

Key words: Low energy electron diffraction (LEED), Low energy electron 
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1 Introduction 



Low-energy electron diffraction (LEED) has been a workhorse technique for 
modern surface science[l]. Experimental measurement of the intensity of the 
diffracted spots as a function of the electron's accelerating voltage (energy), I- 
V analysis, combined with full dynamical calculations of the electron scattering 
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has become over the last two decades one of the most successful techniques for 
surface-structure determination. Although the development of synchrotron- 
radiation-related techniques, recent advances in scanning-probe microscopies, 
and the advent of efficient ab-initio formalisms have decreased its applica- 
tion, LEED remains a very powerful tool. The mature technique routinely 
determines the surface geometry with an accuracy well below a tenth of an 
Angstrom with error limits that are well understood [2]. 

Ruthenium has attracted interest for its catalytic properties [3]. The clean hep 
Ru(0001) surface does not reconstruct. However, like most unreconstructed 
metal surfaces, the topmost layer of atoms relaxes inward towards the sec- 
ond layer. Both early LEED studies [4] and more recent ones [5] reported that 
the first interplanar spacing, di, was contracted by 2% relative to the bulk. 
However, ab-initio calculations by Feibelman[5,6] and Xu et al.[7] found a 
much larger relaxation, between 3.5-4%. This experiment/theory disagree- 
ment prompted suggestions that the smaller experimental relaxation may re- 
sult from hydrogen adsorption [5], although this possibility was disputed by 
others [8]. In a subsequent surface x-ray diffraction (SXRD) study, Baddorf 
et al.[9] found a somewhat larger interplanar spacing and only small changes 
in d\ upon H adsorption. Given this history, the Ru(0001) surface serves as 
a model system to test the accuracy of both experimental techniques and 
ab-initio calculations. 

As Fig. 1 illustrates, the perfect (step-free) Ru(0001) surface exhibits three- 
fold (p3m) symmetry. Thus, diffraction from a single terrace should give a 
3-fold pattern. However, most diffraction techniques show six-fold symmetry 
(p6) in the corresponding diffraction pattern[10,ll]. This apparent change 
in symmetry results from the different stacking on neighbor terraces, with 
termination ...BABA and ...ABAB respectively (where A and B indicate the 
two possible stacking positions of the hep crystal). As Fig. 1 also shows the 
surface terminations of adjacent terraces are rotated by 180° with respect to 
each other; therefore, diffraction from adjacent terraces yield 3-fold patterns 
that are rotated 180° with respect to each other. Thus, averaging over several 
terraces results in the commonly-observed six-fold symmetry. 

In this work we revisit this surface to address the discrepancy between the- 
ory and experiment described above. We use the fine spatial resolution of 
low-energy electron microscopy (LEEM) [12,13,14] to obtain a new set of I- 
V curves from a single terrace on the Ru(0001) substrate. In contrast, in 
traditional LEED measurements, the electron beam illuminates large areas 
comprising many terraces. Obviously, avoiding symmetry averaging should 
enhance the reliability of the I-V analysis. Indeed, our best-fit structure gives 
a first interplanar spacing, d\, that is in excellent agreement with theoretical 
predictions within the carefully determined experimental errors. 
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Fig. 1. (a) Schematic illustration of two adjacent terraces on the Ru hep surface 
separated by an atomic step. The arrows point from a surface atom to its three 
closest neighbors in the layer below. Across adjacent terraces, the arrows are ro- 
tated by 180°. (b)-(c) LEED patterns at 32 eV acquired from two adjacent terraces 
showing the three-fold symmetry of the surface and the rotation upon crossing a 
monoatomic step. The observed spots correspond to the specular and the integer 
(10) beams. 

2 Experimental and Calculations Details 

The experiments were done in a commercial LEEM system with a base pres- 
sure of 7xl0~ n torr. The instrument has a lateral spatial resolution of 8 nm 
and uses a heated LaB 6 crystal as an electron source. The Ru(0001) substrate 
was cleaned by repeated cycles of oxygen adsorption and flash-annealing to 
1700 °C. The crystal contains a large area, about 200 ^m wide, with a low 
density of atomic steps. In this region, terraces more than 5 /im wide can be 
found routinely. The video-rate acquisition system consists of a Peltier-cooled 
CCD camera whose output is digitized by a commercial digital video system 
at 30 frames of 720x480 8-bit pixels per second. The Ru(0001) surface was 
aligned perpendicular to the electron beam using two independent tilt axes of 
a precision manipulator. The substrate/electron-beam alignment was further 
optimized by adjusting the microscope lenses so that the equivalent diffrac- 
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tion beams appeared/disappeared at exactly the same electron energy (i.e., 
the Ewald sphere was centered). 

A major advantage of LEEM is the ability to image a surface and then obtain 
diffraction information from selected and well-characterized areas[12,13,14]. 
This capability is illustrated in Fig. 2. Two curved monoatomic steps, which 
appear dark because of a phase-contrast mechanism [12, 15], are observed within 
the field of view of the left image. An area about two microns wide is selectively 
illuminated by inserting a smaller aperture into the illumination beam (right 
image). By changing the power of the imaging lenses, the diffraction pattern 
from this small area is obtained (Fig. 1). Since this region contains only a 
single terrace (no steps), the LEED pattern is three-fold symmetric (Fig. lb). 
The LEED pattern obtained from a region on an adjacent terrace is rotated by 
180° (Fig. lc), consistent with the change in stacking between adjacent layers 
in an hep crystal. Another way to show the symmetry change when crossing 
substrate steps is to use an aperture [12] to selectively form an image from a 
non-specular diffraction spot [16]. In this dark-field imaging mode, only areas 
that diffract electrons into the selected spot are imaged bright in the real-space 
image. Thus, in an image (see Fig. 3) formed from an integer diffraction spot, 
adjacent terraces are alternately bright and dark. This dark-field contrast is 
similar to the one reported for wiirtize (0001) surfaces[17]. 

For diffraction measurements, a LEEM apparatus has several advantages over 
a conventional electron diffractometer. The electron path is well shielded 
from external magnetic fields, making very low-energy measurements possi- 
ble (VLEED[18]). Since there is no direct optical line of sight from sample 
to the detection screen, high-temperature measurements are straightforward. 
The LEEM geometry, where incoming and outgoing electrons are separated 
by a magnetic prism, allows measurement of the specularly reflected beam 
(i.e., the (00) beam) at normal incidence. The magnification of the LEED 
pattern can also be adjusted. Furthermore the spots do not move when the 
beam energy is changed. Unlike in conventional diffractometers, most current 
LEEM instruments do not perform energy filtering. Thus, secondary electrons 
are imaged together with the elastically scattered electrons. However, the sec- 
ondary electrons can be easily removed by image processing because their 
contribution is spread out over a large area in reciprocal space, as described 
below. 

The intensities of the (00), (01), (10) and (11) beams were obtained simulta- 
neously by recording a sequence of images while changing the electron energy. 
The crystal was maintained at 456 K. To use the full dynamic range of the de- 
tection system and maximize the signal-to-noise ratio, the beam-energy range 
of 1 to 353 eV was split into three ranges (1-70 eV, 60-133 eV and 115-353 eV). 
Within each range, the electron current was adjusted to nearly saturate the 
image collection system at the most intense conditions. The different sections 
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Fig. 2. (a) LEEM image of the Ru surface. The field of view is 15 fim, and the 
electron energy is 2.5 eV. The black dashed lines mark the two atomic steps present 
in the field of view (detected as faint gray lines), (b) LEEM image of the same area 
where a small aperture limits the electron beam to a 2.1 fim diameter region on a 
single terrace. A diffraction pattern is formed from this selected area by changing 
the power of the microscope lenses. 



of the I-V curves were then rescaled with the same scaling factor for all the 
beams, by matching the intensity of overlapping scan regions. A plane was fit 
to the intensity of a box around the diffraction spot. This plane, which con- 
tains the secondary electron contribution, was then subtracted from the box. 
This procedure gives the same results as subtracting the average intensity 
along the box's perimeter, except where secondary electrons impose a large 
gradient on the background. The three (10) and (01) beams, as well as the 
six (11) beams, which are symmetry equivalent [19], where averaged together. 
Given that the (10) and (01) beams are not equivalent, we measured a total 
energy range of 1109 eV. The experimental I-V curves are shown in Fig. 4. 

Full dynamical I-V calculations were performed with a modified version of 
the Van Hove-Tong package[l,2,20,21]. The surface was modeled by stacking 
two Ru atomic-planes on top of a Ru(0001) bulk. Well-converged values for 
both the number of beams and the number of phase shifts (l max — 8) were 
employed. In the structural search only normal displacements for the first two 
Ru atomic planes were considered (di and ofo)- The resulting parameter space 
was fully explored by calculating the I-V curves over a fine 2D grid where 
d\ and d 2 were swept over wide ranges. The experiment-theory agreement 
was quantified via Pendry's R-factor [22] (Rp), while the error bars on each 
parameter were obtained from its variance: AR P = R P)min x yJ(8Vi/AE), 
where V* gives the optical potential and AE corresponds to the total energy 
range analyzed. Correlations between both parameters were taken into account 
for the error-limits estimation. Non-structural parameters such as the muffin- 
tin radius (r m ), the optical potential (Vj) and the Debye temperatures at the 
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Fig. 3. Dark-field image of the Ru(0001) surface. The field of view is 15 //m. The 
electron energy is 40 eV. An aperture blocks all the diffracted electrons except those 
corresponding to one of the first-order diffraction spots: only areas that diffract elec- 
trons into the selected spot are imaged bright. Thus, the image intensity alternates 
between bright and dark when monoatomic steps are crossed. 

surface planes were also optimized (©#, ©d an d Q^ lk for the last layer, the 
next-to-last layer and the bulk, respectively), although they had no impact on 
the final surface geometry. The best-fit parameters are shown in Table 1. 



3 Results and Discussion 

The experimental I-V curves are shown in Fig. 4a. As is obvious from the data, 
the (01) and (10) beams are not equivalent. The best fit to the experimental 
data, shown on the same figure, was found for d x = 2.065 ± 0.02 A and d 2 = 
2.14±0.025 A. With respect to the room temperature bulk value, <4 = 2.141 A, 
these interplanar spacings give relaxations of 5d x = —3.5(0.9)% and 5d 2 = 
0.0(1)% (where the number in parenthesis indicates the error in percentage 
of the bulk value). The dependence of R P on d\ is presented in Fig. 4b. The 
agreement between the experimental and simulated I-V curves is excellent, as 
reflected by an Rp minimum value of Rp ym in — 0.077. The interplanar spacing 
of the topmost two layers for the optimized structure is in excellent agreement 
with the ab-initio calculations. That is, we find that 5d\ = —3.5(0.9)% while 
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Fig. 4. (a) Comparison between experimental (solid lines) and calculated (dashed 
lines) I-V curves for the Ru(0001) surface. The total incident-beam energies spans 
a range of 1109 eV. (b) Rp vs. first interlayer spacing, d\. For each fixed d\ value 
in the graph, Rp has been minimized as a function of c#2, thus taking correlation 
effects into account. Gray horizontal line gives the R-factor variance, ARp, from 
which the error bars are estimated. 

theory finds 3.5-4% [5,7]. However, our structure differs somewhat from the 
previous experimental structures, which found 5d\ = —2% [4,5,8] and 5d\ = 



Given the complexity of the LEEM optics, misalignment of the incidence angle, 
particularly changes with energy, would introduce an additional error in the 
derived structure. We note that the intensities of the individual equivalent 
beams were well matched over the full energy range. This observation strongly 
suggests that the incidence angle did not change with energy. In addition, 
we explored the possibility of a systematic misalignment by optimizing the 
incidence angle in the R-factor analysis. We found that R p had a pronounced 
minimum only at normal incidence and not for any off-normal condition. The 
error bars of this analysis, even with the additional fitting parameter, was 
essentially the same as the best-fit structure of Fig. 4 and Table 1. Thus, we 
conclude that the error in the best-fit structure due to the incidence angle is 
negligible. 

We note that we report results at 456 K. We investigated potential effects of 
this somewhat elevated temperature in two ways. First, we acquired a com- 
plete, independent I-V data set at room temperature. Analysis of this data 
produced the same structure as the best-fit structure at 456 K (Fig. 4 and Ta- 



2.4(0.4)%[9]. 
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ble 1). Second, we optimized the bulk-Ru lattice constant in the R-factor anal- 
ysis of the 456 K data. The minimum R p occurred at the room-temperature 
Ru lattice constant (within ±0.2%) and, again, the error bars were nearly the 
same as the analysis without this additional parameter. Thus, we conclude 
that the first interlayer spacing of Ru is contracted by 3.5% both at room 
temperature and at 456 K. 

We next discuss the advantages of determining a surface structure from a single 
terrace instead of using multiple terraces. First, for single-terrace data a larger 
energy range, AE, can be obtained since symmetry-inequivalent beams are not 
averaged together. This larger energy range improves the reliability of the fit. 
Indeed, our energy range is larger than that used in previous LEED analysis[5]. 
To evaluate this potential error source, we performed an independent R-factor 
analysis of our experimental data where the experimental (10) and (01) beams 
were averaged. Compared to the single-terrace data set, the best fit for this 
symmetrized data set shows a slightly smaller interplanar contraction but a 
significantly larger error bar, 5d\ = —3.3(1.4)%. With this larger error, this 
result would be compatible with the previous experimental results. 

An additional uncertainty in multiple-terrace data is that the two terrace 
types are usually taken as being present in equal abundance. But this is not 
always the case. Ru(0001) has a marked tendency to form "double-steps," 
that is, one terrace type is considerably wider than the other terrace type. 
The narrow terrace type is bounded by two closely spaced ("double") steps. 
Which terrace type is more abundant depends upon the local orientation of 
the steps. Even for an averaging technique, such as obtained from a traditional 
diffractometer, the two terraces types may not be in equal abundance, giving 
an extra parameter to fit and increasing the error bars. This effect might 
explain the experimental observation of three-fold patterns on Ru surfaces by 
X-ray photoelectron diffraction [23]. In fact, the change in symmetry has been 
employed to monitor the presence of double steps[10]. Double-steps may also 
introduce a non-negligible correction to the usual assumption of incoherent 
beam mixing in the theoretical I-V simulations, further reducing the reliability 
of the R-factor analysis. 

The final point concerns the role of adsorbed hydrogen in the experimental 
structures 1 . We note that our structure, which was obtained from a crystal 
above room temperature, agrees extremely well with ab-initio calculations. 
Since there is no disagreement, we find no need to invoke impurity effects. 
To our knowledge, this is the only experimental result which corroborates the 



The LEED analysis of Menzel et al. as well as the SXRD analysis of Baddorf et al. 
concluded that H adsorption was not responsible for the experimentally determined 
Ru(0001) first interlayer spacing[8,9]. However, the calculations of Feibelman et 
al. [5] show that H adsorption significantly decreases the surface contraction. 
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Structural Parameters [A] 


Non-structural Parameters 


di 2.065 ± 0.02 

o 1 /in _i_ n c\ok 
(12 Z.14U ± U.Uzo 

dtuik 2.141 


r m 2.00 bohr 
Vj -z.d tj ' [evj 
9^ 350 K 
@ 2 D 400 K 
e^' fc 600 K 



Table 1 

Optimized structural and non-structural parameters deduced in this work. See text 
for further explanations. 

3.5% contraction deduced from the hydrogen- free DFT calculations. 



4 Summary 

Using a LEEM instrument we have performed, for the first time, a LEED I-V 
analysis from a single atomic Ru(0001) terrace. The experimental I-V curves 
were excellently fit by full dynamical calculations. In the best-fit structure, 
the topmost Ru layer is relaxed inward by about 3.5%, in excellent agreement 
with ab-initio calculations. We suggest that determining structures from small 
areas has clear advantages even for surfaces that consist of a single structure 
(phase), such as Ru. 
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